Quantum Dynamics of Spins Coupled by Electrons in ID Channel 
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We develop a unified theoretical description of the induced interaction and quantum noise in a 
system of two spins (qubits) coupled via a quasi-one-dimensional electron gas in the Luttinger liquid 
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Recently, there has been much interest in coherent 
quantum dynamics of coupled two-level systems (qubits) 
for quantum information processing. Realizations are 
sought such that qubit-qubit interactions can be exter- 
nally controlled over short time scales of quantum "gate 
functions," in the parameter regime ensuring that relax- 
ation and decoherence are negligibly small over a large 
number of gate cycles. There have been several pro- 
posals for qubit systems in semiconductor heterostruc- 
tures, with direct coupling [1], typically via shared elec- 
tron wave functions, or indirect coupling, specifically via 
excitations of the conduction electron gas [2] . In the lat- 
ter approaches, the medium that induces the indirect in- 
teraction, can also act as a "heat bath" resulting in re- 
laxation and decoherence. Usually, strongly correlated, 
low-temperature conditions have been assumed [3] in or- 
der to ensure high degree of coherence. In this work we, 
for the first time, develop a unified theoretical derivation 
of the Ruderman-Kittel-Kasuya-Yosida (RKKY) type [4] 
induced interaction incorporating the description of re- 
laxation effects resulting from the electron gas "bath." 

A recent experiment [5] on coupled quantum dots has 
demonstrated the realizability of indirect interaction for 
control of two-qubit dynamics. Several experimental se- 
tups [6] suggest that a quasi one-dimensional (ID) chan- 
nel geometry for the conduction electron gas is promising 
for quantum measurement required for quantum comput- 
ing. Furthermore, there is experimental evidence [7] of 
Luttinger liquid behavior in electron transport in quasi 
ID structures. Therefore, we are going to consider the 
ID-channel qubit-qubit coupling via indirect RKKY in- 
teraction mediated by Luttinger liquid of electrons, and 
we assume spin-1/2 qubits. 

It has been commonly accepted [8] that the ground 
state excitations of a ID interacting electron gas within 
the Luttinger liquid model can be described by the fol- 
lowing Hamiltonian, 



(1) 



The phase fields, 4> c (s)( x ) an d 9 c (s){ x )i with sub- 
script indices describing charge and spin degrees 
of freedom, respectively, obey commutation relations 



[d x <f>c(8)( x )iQc(s){ x ')] = 2mS(x— x'). Consequently, 9 and 
d x 4> can be viewed as canonical variables. Here and in 
what follows we set Ti — 1 and fee = 1. The Hamiltonian 
has a simple additive structure as a result of spin-charge 
separation in ID systems, with the charge and spin den- 
sity waves of the liquid having, generally speaking, differ- 
ent velocities, v c = v F /g c and v s = v F /g s , respectively, 
where v F is the Fermi velocity. The constant g c > ac- 
counts for the electron-electron interaction and is related 
to the parameters of the Hubbard model [9] as follows: 
g c ~ (1 + U/2E F y 1 / 2 , where E F = v F k F /2 is the Fermi 
energy, k F is the Fermi momentum, and U is the effective 
interaction between the electrons, U ~ e 2 /a, where a is 
the short distance cutoff, a ~ kp l . Also, we assume rota- 
tional symmetry, SU(2), in the spin space [9-11], which 
implies that g s = 1. 

The localized magnetic moments (spins) are coupled 
to conduction electrons via the contact interaction, 



Hint = 



3 



(2) 



Here j = 1,2 labels impurity spins Sj positioned at Xj, 
Jj are the exchange coupling constants, and s(x) is the 
local electron spin density. The spin density can be ex- 
plicitly expressed in terms of the Luttinger phase fields, 
see [12,13], 



—!—!L 4- — cos(2k F x + 6 C ) cos S , 
2n ira 



(3) 



tt a 



[Huy cos 8 S + a x cos{2k F x + 8 C )\ , (4) 



where <J x . VtZ are the Pauli matrixes. The Luttinger liquid 
description of the problem is generally valid in the "hy- 
drodynamic limit" of spin separations x — \x\-x-z\ a. 

Our goal is to obtain an effective description of the 
dynamics of the system of two spins, with electronic 
degrees of freedom integrated out. We first consider 
the equilibrium partition function of the system, de- 
fined as Z = Tr[exp(-/3if)], where = l/T is the 
inverse temperature and H = H + H- mt . The parti- 
tion function, Z , can be expressed in terms of the spin- 
and temperature-dependent effective action 5 c ff, see [14], 
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Z = J DS1DS2 exp [— <S e ff(Si, S2)]. As usual, the evalua- 
tion of the effective action can only be carried out pertur- 
batively in the couplings Jj of the spins to the electrons. 
The leading non-vanishing contribution is generated by 

(1/2) J Q dT 1 dr 2 (T i?int(n)-tfint(T2))ir , where T stands 
for Matsubara time ordering [15], and the equilibrium 
averaging is taken with respect to the non-interacting 
Hamiltonian Hq. The perturbative approach is, gen- 
erally speaking, invalid at sufficiently large Matsubara 
timescales, i.e., low temperatures, as evident from the 
perturbative RG analysis carried out in [11,12]. At suf- 
ficiently low temperatures, spin dynamics results in a 
nontrivial strong coupling fixed point, i.e., Kondo effect. 
For quantum information processing, we are interested 
in qubits that retain their localized spin-1/2 indentity. 
Therefore, we limit our consideration to temperatures 

T > Txondo- 

The resulting Matsubara action is <S c ff = iSecrry +<?seif • 
Here iSecrry is the Berry action term for the spins. It will 
be addressed later. Presently, we focus on the self-action 
for the spins, given by 



Airvp 



+ 



C2{g) cos (2kFx) 



-v- 



(7) 



where C x = {2n)-a- 1 g9 j dz(l + z 2 )- 1 / 2 (g 2 + z 2 )-9/ 2 and 
C 2 = 4(2tt)- V _1 r(l - g ) sin [(tt/2)(1 - g)}. The first 
term in (7) corresponds to interaction between the spins. 
The interaction is oscillatory and decays as a power law 
x~ 9 . This result is consistent with [13]. 

The second, w„-dependcnt term in (7) corresponds to 
relaxation of the spins. To demonstrate this property, 
we perform a transition to the real time dynamics of the 
spins according to analytical continuation rule, see [16]. 
We introduce a standard Keldysh contour with forward 
and return branches time-ordered and anti-time-ordered, 
respectively. In the Keldysh representation, the effective 
action, (5), reads 



~2~ J 2^ 
+(Si 



Sf(w)x(w,a;)Si(-u;) 
S 2 ) . 



S[( W )x(w,0)S 2 (- W ) 
(8) 



J 2 r 

<S se if = |x(^n,0)Si(iw„) • Si(-iw„) + 



x(u n ,x)Si(iu> n ) • S 2 (-iw„) + (Si <-> S 2 ) 



In (5), xi^nix) is the spin-spin correlation function of 
the electron gas (in the imaginary time rcprescnation), 
and oj n — 2im/ [3 are the Matsubara frequencies. For 
simplicity, in (5) and in the following we assume that 
Ji = J 2 = J. The correlation function x(w„,a;), which 
generally speaking is a tensor, (si(— iu n , x)sj(iu) n , 0)}, in 
the SU(2) symmetric case reduces to a scalar function, 
X{u n ,x) = {s z (—iuj n ,x)s z (iuj n ,0)). It can be evaluated 
by using (1) and (4). Assuming that the electron gas is 
dense to satisfy Ef ^> one obtains (with g = g c ), 



Here Sj=i ; 2 are two-element column vectors, with spin- 
vector operators as elements, such that their transposes 
areSf = (Sj c ,S, 9 ). The "classical," S, c , and "quantum," 
(5) Si q , components of the spin Si are calculated as Fourier 
transforms of the following combinations of spin opera- 
tors on the forward (/) and return (r) branches of the 
Keldysh contour, S ic = (Sj^ + Sj r )/2 and S q = Sif — S ir . 
The response (correlation) function x( w > x ) is then a 2 x 2 
matrix, which can be expressed in terms of the retarded 
and advanced response functions, x R an d X A L > 



4irv 2 F 



exp(- 



') + 



cos (2kFx) 

(2 7 r)9+ 1 a 1 -9 J (x 2 + v 2 f t 2 ) 1 / 2 (x 2 +v 2 F T 2 /g 2 )9/ 2 



V F 

dTexp (— iui n r) 



(6) 



The correlation function in (6) contains two contribu- 
tions, one due to forward scattering, peaked at wavevec- 
tor (7 = 0, and another due to backwards scattering, 
peaked at q = ±2fci? . As expected, the forward scat- 
tering contribution is g-independent. 

In the perturbative regime considered, the dynamics 
of the spins is slow and controlled by small parame- 
ter J 2 /v 2 F . Therefore, we can use the small- frequency 
asymptotic form of the correlation function x(w„, x) for 
\u n x\/v F < 1, 



C\(g) cos (2ftFx) 
VFa x ~ 9 x9 



X 



X A 
x R x K 



(9) 



The retarded and advanced response functions are re- 
lated to the Matsubara response function via the ana- 
lytic continuation iuu — > u ± iS. In thermal equilibrium 
X K can be expressed in terms of the retarded and ad- 
vanced components by using the fluctuation-dissipation 
theorem, x K = coth (/3w/2)(x fl - X A )- 

Let us consider first the noninteracting case, g = 1, 
and later we will extend the results to g ^ 1 . For nonin- 
teracting electrons, the response function corresponds to 
an Ohmic heat bath, with Ci(l) = 1/(4tt) and C 2 (l) = 1 
in (7). Upon Fourier transform, (8) yields several terms. 
Those containing products Si c ■ Sj q represent interaction 
between spins, while the Si C ■ Sj q and Si q ■ Sj q terms 
are responsible for energy dissipation and pure dephas- 
ing (decoherence), respectively. The dissipative (time- 
derivative) terms are small and can be neglected here. In- 
deed, since S, ~ 0(J 2 /v F ), the S-dcpcndent terms in the 
action in (8) are of order J 4 jv F . The resulting dynamics 
of the spins is governed by the action S f 
where 



-'self 



°int ' °dec' 



2 



J cS J dt[S lc (t)-S 2q (t) + S lq (t)-S 2c (t)} , (10) 



S ^ cc = 2~7r J dtldt2 J dcjuJC0S M f 2 - h)\ x 
coth {^^j {2Si g (ti) • Si,(t 2 ) - [1 + cos(2fc F x)] x 
Si,(ti) • S 29 (i 2 ) + [Si, <- S 2 ,]}, (11) 



with 



2 cos (2fc F a;) 

Jeff = J Z 

Attvfx 



and 7 = 



J 2 



2ttv 2 f 



(12) 



Equation (10) corresponds to coherent interaction of the 
spins according to the scalar Si • S 2 coupling. Equation 
(11) represents quantum noise resulting from thermal 
fluctuations of the Luttinger liquid. The noise, generally 
speaking, is colored. Moreover, the noises experienced by 
the two spins are correlated. 

The action in (11) can be simplified further if we recall 
that the dynamics of the spins is slow, i.e., S iq (t 2 ) = 
S iq {ti) + 0(J 2 /v 2 F ). Then, the self-action can be replaced 
by an instantaneous action, by setting Sj 9 (ti) = Sj 9 (f 2 ) 
in (11), which now corresponds to the white noise source. 
The full action of the spin system can now be written as 



S K =Sg my -J dtH K (Si.Sa), 



(13) 



where the Berry action plays the role of a velocity x mo- 
mentum term in the Largangian-Hamiltonian transfor- 
mation [17]. The "Keldysh Hamiltonian" in (13) is 



H K = — J e g (Si/ • S 2 / — Si r • S 2r ) 



(sin 4> S lq - cos cf> S 2g ) 2 + (Sig <-> S 2g ) 



(14) 



where < (f> < ir/2 is defined by sin (2</>) = cos 2 (fc F x). 

The density matrix of the two-spin system evolves ac- 
cording to the "Schrodinger equation" p — —iTk{H k p) 1 
where Keldysh time-ordering implies that the "forward" 
operators, with subscript /, are positioned to the left of 
the density matrix, while the "return" operators, labelled 
by r, are to the right of the density matrix, 

p = iJ e s[S 1 -S 2 ,p]-4ryTx 

( [ sin 0Si - cos <j>S% , [ sin <j>S? - cos <t>S% , p] ] 



a—x,y.z 

+ (Si ^S 2 ) 



(15) 



This expression contains both the mediated interaction 
and dephasing due to the electron environment. 

It is instructive to illustrate the dynamics of the two- 
spin system, described by (15), for a particular initial 
state of the system, ||J.). Without the quantum noise 



term, proportional to 7T in (15), the Si • S 2 interaction 
would split the singlet and triplet spin states. As a re- 
sult, the system would oscillate between the |t|) and ||t) 
states with frequency determined by the singlet-triplet 
energy gap. The effects of the noise include damping of 
these oscillations, as illustrated in Fig. 1. Furthermore, 
for t > 0, the system subject to the noise will no longer 
remain in a pure quantum state. The departure of the 
resulting mixed state from a pure state can be measured 
by deviation of Tr[p 2 (i)] from the pure-state value of 1, 
as illustrated in Fig. 2. We point out that effective evo- 
lution equations that involve only commutators, linear 
in the density matrix, on the RHS, typically fail to re- 
produce thermal equilibrium at large times. Instead, as 
seen in Figs. 1 and 2 — note the assymptotic values — 
the fully random mixed state is obtained. Therefore, the 
present approximation should not be used beyond the re- 
laxation time defined by (15), namely, it only applies for 
t < l/[T(J 2 /v F )], and the theory is therefore applicable 
in the regime of interest for quantum computing applica- 
tions, for short and intermediate times, because both fac- 
tors in the denominator, T and J 2 /v F , are small. Similar 
results have been obtained in Ref. [18], studying medi- 
ated interaction and decoherence due to noninteracting 
electron gas. 

Finally, we extend our results to the interacting case, 
g 7^ 1. Modification of the interaction term, (10), is 
straightforward and comes from the first term on the 
RHS of (7). Equation (10) applies with 



Jeff(3) = Ci( 5 )J ; 



2 cos(2fc F a;) 



(16) 



Modification of the decoherence term, (11), is less ob- 
vious. Consider for simplicity a single spin situation, 
when the first term on the RHS of (7) can be ommitcd, 
and in the second term we can put x = 0. The \u> n \ 
term then in the brackets in (7) produces the same, g- 
independent, contribution to the decoherence rate as in 
the noninteracting case. The \to n \ 9 term in the brackets 
requires careful consideration. Indeed, we are primarily 
interested in the repulsive, U > 0, Hubbard-model inter- 
action, i.e., < g < 1. The two-spin contribution them 
yields a divergence: After setting Sj 9 (ii) = S iq (t 2 ) in 
(11) and integrating over t\ — t 2 , the frequency integral 
/ duju> 9 S(uj) coth (Poj/2) is divergent for g < 1. The ori- 
gin of this divergence is actually not in the instantaneous 
assumption but instead it can be traced back to too lib- 
eral a use of the large-/? approximation in our evaluation 
of the response function in (6), and specifically, extend- 
ing the limits of the r integration from -co to 00, while 
we should have integrated from —/3 to (3. One can show 
that 8{u) should be regularized as, e.g., sin 2 (/3w)/(7r/3w 2 ), 
which eliminates the unphysical small-frequency diver- 
gence. As a result one obtains the dephasing rate 



3 



V v g F a 1 9 J 

where C3 <~ 1 for g ~ 1. Equation (17) represents mod- 
ification of the Korringa law for spin relaxation in the 
interacting ID electron gas. A precise determination of 
the coefficient C3 is out of the scope of the present work 
and will be dealt with in a subsequent analysis based on 
a fully non-equilibrium description of the system. 

In summary, our main result, (15), is the first the- 
oretically derived dynamical equation that incorporates 
both the coherent RKKY-type induced interaction and 
the effects of quantum noise due to interacting electrons 
in ID conduction channel. The approximations and as- 
sumptions involved, limit our results to temperatures in 
the range Tk < T <C Tp. Furthermore, to have the 
noise term small, the temperature should be actually in 
Tk < T < vp/x, where x is the qubit (spin) separation: 
restoring the constants earlier set to 1, the upper bound 
here is Uvp/ksx. The relative strength of the interac- 
tion vs. noise terms can be also controlled by positioning 
of the qubits, owing to the oscillatory dependence on x, 
typical for RKKY-coupled systems. 
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FIG. 1. Lower curve: The probablity, P|t|),|u>' to find the 
two-spin system in the |f J.) state as a function of time, for a 
convenient set of parameter values, J e g = 1, 7T — 0.0125, 
k F x — ir. Upper curve: The sum of the probabilities to 
find the two-spin system in the states |||) or |||), given by 
P|Ti),ITI> + PUT), lit)- The large-time limiting values of these 
probabilities are 1/4 and 1/2, respectively. 




FIG. 2. Solid line: The quantity Tr(p 2 ), with the initial 
value of 1 corresponding to a pure state, as a function of 
time for the same parameter values as in Fig. 1, with posi- 
tive J e f{ — 1. Dashed line: Parameter values illustrating the 
case of negative J cf f — —2, with 7T = 0.0125, k F x = tv/2. 
Dotted line: The deviation form a pure state (the quantum 
noise effect) is present even when the leading-order induced 
interaction is J e s — 0, with 7T = 0.0125, k F x = n/4. The 
curves for J" e fi / have weak oscillations superimposed on 
the decay. The large-time limiting values are 1/4. 
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